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If a primordial magnetic field (PMF) is present during photon decoupling and afterward, a finite 
neutrino mass can affect all modes of the CMB. In this work, we expand on earlier studies of the 
scalar mode effects by constructing the vector and tensor mode equations in the presence of massive 
neutrinos and a PMF. We compute the power spectrum of the various modes in an illustrative 
example and find that the neutrino mass can significantly affect the vector and tensor modes when 
a PMF exists, while the effects are negligible for no PMF. The most prominent result of the present 
analysis is the behavior of the EE component of the tensor mode at low multipoles. For massive 
neutrinos the EE mode can become comparable to the observed primary anisotropy. Therefore, if 
and when the EE mode power spectrum is measured at low multipoles the possibility exists to place 
a strong constraint on the sum of the neutrino masses. 

PACS numbers: 98.80-k, 98.80.Jk, 97.70.Vc, 14.60.Pq, 52.27.Ny, 98.65.Dx, 98.70.Vc 



I. INTRODUCTION 



The possibility that a primordial magnetic field (PMF) 
could affect the CMB power spectrum is a subject of con- 
siderable recent interest, e.g. P, H, H, Elj Hi • Such effects 
arise in two ways. Baryons are affected by the Lorentz 
force from a PMF. They then influence the CMB spec- 
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trum indirectly through Thomson scattering. Moreover, 
the Lorentz force from a PMF also has a significant ef- 
fect on the later development of large-scale structure. 
Among other things a PMF could explain the possible 
power excess observed by ACBAR Q for large multi- 
poles {I > 2500) in the CMB power spectrum. It could 
also explain the BB mode CMB anisotropies observed by 
CBI [S] as well as the origin of the magnetic field inferred 
from the observed [9l.[l0l.flll.[T^ polarization of light from 
galactic clusters. Such studies suggest that a PMF may 
indeed exist, and that it is worthwhile to examine other 
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possible consequences of its existence. In previous work 
|13| . we deduced neutrino mass effects from the scalar 
mode of the CMB power spectrum using the WMAP-3yr 
data when the PMF effect is also taken into con- 
sideration. Indeed, the neutrino mass constraint from 
the scalar mode of the CMB is the recent focus in the 
literature [ij, [l^, In this work, we expand on this 
discussion by introducing new studies into the neutrino- 
mass effects on the vector and tensor modes of the CMB 
as well as the scalar mode. 

The CMB power spectrum induced by the PMF sep- 
arates into 3 parts, i.e. the scalar, vector, and tensor 
modes. The scalar mode affects the power spectrum at 
large angular scales I < 10, while the vector mode primar- 
ily affects I > 1000. In previous work [l^, we deduced 
neutrino mass constraints from the scalar mode of the 
CMB power spectrum using the WMAP-3yr data. We 
showed that a PMF can decrease the upper limit on the 
neutrino mass. Our ultimate goal, however, is to con- 
strain the neutrino mass using all PMF modes and using 
all available observations including all multipoles. How- 
ever, before calculating the neutrino mass constraint, we 
need to derive a framework in which to study the effects 
of the neutrino mass on the CMB power spectrum in the 
presence of a PMF. Hence, the purpose of this work, is to 
construct the relevant equations for the effects of massive 
neutrinos on the vector, and tensor modes in addition to 
the scalar modes of the CMB power spectrum. With 
this framework in place, we then study the possible neu- 
trino mass effects in the context of an illustrative model. 
We show that the neutrino mass can cause the EE mode 
at low multipole to become comparable to the observed 
primary anisotropy. We also find that the BB mode is af- 
fected by the neutrino mass and may be detectable if the 
gravitational wave contribution is not too large. Thus, if 
and when the EE and BB mode power spectra are mea- 
sured the possibility exists to better constrain the sum 
of the neutrino masses. 



A. Boltzmann equations and anisotropies 

To begin, we choose metric perturbations in terms of 
the scale factor a and the conformal Newtonian-gauge 
metric parameters ^ and <& of the form. 



(5,900 = -2a2*g(°), Sg,j = 2a^'i>Q'-°^S.. 



for the scalar mode, and 



(1) 



(2) 



for the vector mode, while the tenor mode metric pertur- 
bation is given by 



(3) 



Here, Q^°\ Q^f and Qf are the harmonic modes for 
the scalar, vector and tensor components, respectively, 
while the H^^ , h!^^ are vector and tensor metric pertur- 
bations, respectively. Here and throughout we use the 
same notation as in Ref. J^]. Our gauge for the vector 
mode, however, is different from that of Ref. [l9|, and 
the same as in Refs. [1, i-e- the perturbations for the 
vector mode are in space-space components rather than 
space-time components. This is because it is difficult to 
treat massive neutrinos in the gauge of Ref. [iB] . 

Next, we define coordinates by setting k = 63, and 
defining basis vectors: 



3(±) 



\/2 



(ei±ie2), e. 



(±±) 



;(±) 



3(±) 
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We also write the distribution functions for massive neu- 
trinos as 

f{q,n\k\r^) ^ f{q) + 6f{q,n\k\7j) = fil + e^) , (5) 



II. EQUATIONS 



Here, we derive equations for the perturbations on the 
CMB anisotropies from the various modes in the pres- 
ence of a PMF. Usually, only the scalar mode CMB 
anisotropies have been calculated as in Ref. [17]. How- 
ever, in the presence of a PMF the vector and tensor 
modes @, are also important. Our equations are based 
upon the formulation given in Refs. J^, JJ^]. The equa- 
tions in Ref. [l^ are almost the same as in Ref. [l^ ex- 
cept for the PMF terms. However, a finite neutrino mass 
was not considered in Ref. Hence, in this paper we 
extend the derivation in Ref. to deduce equations 
which include the effects of a finite neutrino mass. 



where Sf is the perturbed part, q is the comoving mo- 
mentum, is its direction, 77 is the conformal time, and 
Qh is the normalized perturbation, with the subscript h 
denoting hot dark matter. Hereafter, we attach subscript 
h to distinguish any quantities Xh which are associated 
with hot dark matter (i.e. massive neutrinos). We will 
confine the subscript v to denote for massless neutrinos, 
i.e. X^. 

Neutrinos do not interact with other particles. Hence, 
they obey the coUisionless Boltzmann equation. We use 
the linearized coUisionless Boltzmann equation to calcu- 
late the evolution of the perturbation, i.e. 



6f = -'^^!^Sf + {q^ + iken%^)^ 
e oq 

'"'^Y,^Sf+iq^ + ^ke^YM% ,(6) 
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for the scalar mode, while the vector mode equation is 



stress TTfi as follows: 



e ■' oq 



and the tensor mode equation is 



5f 



'-^y^^I-^iH^^'^Yfl , (8) 



where e = y^T^ + mJc?, in the same notation as in 
Ref. flTp. To get hierarchial equations similar to Eq. (60) 
of Ref. [19|] , we next expand the perturbation in spherical 
harmonics, 



,(0) 



,(1) 



-(1) 
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'hO 
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Ph 
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q^dqqM^l 



'h2 



a / q dq—f(d\ 
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_4 87r/15 /■ 2 , 9^ jp,(2) 

a / q dq—fQl 

Ph J e 



'h2 



(13) 



where ph and are the energy densities and pressures 
from hot dark matter, i.e. massive neutrinos. Equations 
for the other components are the same as in Ref. [l^ 
except for the gauge difference in the vector mode. This 
gauge difference changes Eq. (61) of Ref. 19] to 



hi 



(9) 



1.771 



= fp(i) - ^4^) 



(14) 



where the olj™"* are amplitudes for the various modes. 
We use the Clebsh-Goldan relation: 



yOy^n 



V(2; + i)(2;-i) 

^{2l + l){2l + i) 



to expand the Boltzmann equations in terms of the Q 



C)(™). 



hi 



A(m) _ qk ( \Jl'^ ~ra^ ^(m) 



^(/ + l)^-m2 ^(^) \ („) 

^7T3 ^'^mj 



(11) 



where, the source terms 5^™^ are given by, 



^In/ ^(0) ^ fee gin/ 
91ng' ''^ g 91ng 



ij^^^ a In/ 
Sing' 



^(2) _ jj{i) d\^! 

0^2 — 



^ 9 In (7 



(12) 



This is very similar to Eq. (60) of Ref. [T^ , except for 
the source terms. This difference is due to the finite 
neutrino mass which makes it impossible to integrate the 
distribution functions easily. The O^^™'' are related to 
the density perturbation (5;i, velocity u/j, and anisotropic 



In order to obtain the power spectrum of the CMB and its 
polarization, we need to expand the temperature pertur- 
bation with spherical harmonics, and expand the po- 
larization fluctuation Q ±iU with spin-2 harmonics [l9|. 

The expansion coefficients are O^™'', and E'^'™' ± iB^^^ 
respectively, and their integral solutions are given as fol- 
lows fii. 



e[°^(r/o,fc) 
21 + 1 



21 + 1 



epWfc) 

21 + 1 

21 + 1 

Bt\rio,k) 
21 + 1 



f + v]/ - $) j 



• (0) -(10) 



drje 



■ (1) -(11) 



(21) 



d7ye--[fF(2)_ij(2)]^-(22) 



Ve / d7jre-^P'-'^h\"'^ 
Jo 



V6 / driTe-^P^"'^/3l"'^ 



(00) 

I 

(15) 

(16) 
(17) 

(18) 

(19) 



where rjQ is the present conformal time, while the ra- 
dial temperature function ™' (x) , the radial E function 



e'f^\x) and the radial B function l3i"''{x) are evaluated 
at X = k{r]o — rf). Here, we have used the anisotropic 

scattering source = (e^"^-V6^2"')/10- The equa- 
tion for the vector anisotropics is different from that of 
Ref. ^] because we are in a different gauge. 



(m) , 
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From Eqs. (|T5)) - (|19p . we can then derive the 
CMB power spectra of temperature and polarization 
anisotropies by constructing correlation function [l^ ; 



{21 + i)2cf 



k^xl"'>im,k)Xr'im,k) , (20) 



(m) , 



where X is Q, E or B. 

The evolution of the perturbation variables is given by 
the Einstein equations [IB] to be: 



(pfSf + P^Ppmf) 



+ S-ipf+pf)vf/k 



fc^(* + $) = -87rGa^(p/4*'^ +Pi^pmf) . (21) 

for the scalar mode, where pf and pf are the energy den- 
sities and pressures from the nonrelativistic (e.g. mat- 
ter) and relativistic (e.g. photon) fluids (following the 
notation in Ref. [T^). ppmf and 7r^|^ are the energy 
density and anisotropic stress from the PMF, which are 
normalized by the photon density p^ and photon pres- 
sure py, respectively, for the following reason. Energy 
momentun tensor for the magnetic field is defined by: 



Stt 



(22) 
(23) 



where we have neglected the electric field because the 
coductance of the early universe is assumed to be very 
large by taking the MHD approximation. Transforming 
Eqs. and (|23p into wavenumber-space and decom- 

posing the traceless part of Eq. ([23]) into scalar, vector 
and tensor mode, we can define the energy density and 
anisotropic stress for the PMF. Since we assume that the 
conductance is infinite, the magnetic field is "frozen in" 
and the time evolution of i? is B cx a^^. This means that 
the energy density and anisotropic stress for the PMF 
grow as oc a~*. To eliminate the effect of expansion of 
the universe, we define the energy density and anisotropic 
stress for the PMF divided p^ (oc a~"') and p^ (oc a^^), 
respectively, as shown in Eq. (|2ip . In our definition, 
therefore, ppmf and ttpmf are comoving variables which 
stay constant of time[l, 0, H, [U. 

The vector mode evolution equation becomes, 

= SnGa^pfirf^ +p^TT%p) 
a •' 



= SirGa 



2 "total 



(24) 



For relativistic components, the pressure Pf is propor- 
tional to a~^. Then, assuming a scaling relation. 



(25) 



and defining the shear as. 



(26) 



we get a simple solution to the vector mode equation 
Eq. jH, 



C2a ^ (Radiation — dominated) 
""^■^ ' C2a~2 (Matter - dominated) . (27) 



Here, the power spectral index a and coefficients ci and 
C2 depend on whether one is treating the radiation or 
matter dominated epoch, as discussed later. 

The shear affects the CMB spectrum through Eqs. 
p6)) and p?)) . If there is no source of anisotropic stress, 
i.e. tt'^P — TTplj^ — 0, Eq. ([27]) has a simple decaying so- 
lution a^^^ — C2a^^ ■ This means that there are no vector 
CMB anisotropies if there is no anisotropic stress because 
the potential decays rapidly. This is the reason why the 
vector mode is usually ignored in CMB theory. However, 
if there is a source of anisotropic stress, a'^^^ decreases 
slowly. 

The tensor mode evolution equation is 



a ■' 



= SnGa 



/2) 
2 "total 



(28) 



Note, that unlike the vector mode evolution [Eq. ([M]) ]. 

there is a linear term, k^H!jp \ on the l.h.s. of Eq. (|28p . 
Assuming that all scales are outside of the horizon, we 
get a solution for the vector mode which is similar to Eq. 
(HZD, i.e. 



C2a (Radiation — dominated) 

-2 



ct(^) = cia""i-^ + C2a"2 (Matter - dominated) . (29) 
This implies that, if there is a primordial magnetic field, 

(2) (2) 

Hrp ' grows even if there is no initial H!^ ' . This passive 
mode was studied in detail in Ref. Q . However, we do not 
consider this mode further here, because we are interested 
only in the neutrino mass effects. 

The essential effects of massive neutrinos have been 
analyzed previously in Refs. [l^ [12, [2^. There are two 
primary effects. First is the free-streaming effect. Out- 
side the horizon, there is no free streaming effect because 
the qk/e term on the r.h.s. of Eq. (fTTj) is negligible. Of 
course, the Os do " free stream" to neighboring hierarchi- 
cal equations once a scale enters the horizon. This free- 
streaming effect , however, decreases as neutrinos become 
non-rclativistic because the qk/e term on the r.h.s. of Eq. 
(fTTj) becomes negligible again. 

The second effect is the change of anisotropic stress 
at the epoch when massive neutrinos become nonrela- 
tivistic. Although anisotropic stress will decrease by the 
effect of neutrino mass, this change is very small without 
a PMF. On the other hand, the existence of the PMF 
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affects drastically the evolution of the anisotropic stress. 
We can understand this clearly from a comparison of the 
initial conditions for massless and massive neutrinos with 
and without a PMF. 



B. Initial conditions for massless neutrinos 

We can derive initial conditions for massless neutrinos 
from the assumption that there are no radiative vorticity 
and anisotropic stress at very early times, i.e. 7rj™^j = 0. 
The result is [2, ji] 



_ _ (1) R-y f 



1 - 



45 {kr]f 



r(2) = 



V 14 4i?^ + 15 



.(1) 



.(2) 



15 {kT]f 



14 4i?, 



15 



15 ' 
(30) 



where P(A;) represents the power spectrum, ttI'"'' is the 
neutrino anisotropic stress normalized by-\/P(fc), R^ = 

P^/pr, Ru = Pvl PR and/9_R = p^+p^- Here cTj^^-* andiJ^^] 
are the initial conditions for the primary mode. These 
are taken to be initially turned off in the PMF mode. 
Hereafter, any quantities denoted with a hat (e.g. X) 
are the variables normalized by the square root of power 
spectrum, P{k). They are very useful variables for the 
illustration of the effects of massive neutrinos to be given 
in Figs. [2] and [3 



C. Initial conditions for neutrinos with finite mass 

The initial conditions for massive neutrinos are a bit 
different from those of massless neutrinos. Expanding 
q/e ~ 1 — rr?a} l(2q^)^ Eq. can be expressed as 



^const. f 3 l"i^a^NF^(m) 



Ph 

^^J-\l-\^,Hl^RmW) 



We here define the effective wave number kef / by 

klff = k^^kl , (32) 
ri 5 \ 4R.. 4- 15 

h — 



1 5 4i?. + 15 



(33) 



conditions for the total anisotropic stress tTj^^L as 



total 



45 {keffj^f 



^PMF 



total 



14 + 15 a' 

(2) 
' total 



Pj_ _ (1) 2fc?7 pi, 
R, a-4 



^(2) 15 {keffV? 



' PMF 



14 ARy + 15 a- 



T,i 



ARy + 15 a- 



(34) 

This condition is also valid for massless neutrinos by set- 
ting = 0, for which fc, 



k. Recall that the a. 



(1) 



and ^ are the initial conditions for the primary mode, 
which should be turned off in the initial PMF mode. This 
equation clearly shows that the neutrino mass effects cor- 
relate with the PMF to leading order in rj^ and would 
affect strongly the PMF mode of the CMB temperature 
and polarization anisotropics. Note, that the expressions 
in Eq. p4|) are only valid for the epoch before the mas- 
sive neutrinos become nonrelativistic, i.e. a < a„r, which 
is given by \15^ 



T 

^ v.nr 

1.5 X 10^3 



(35) 



Thus, a„r = 0.84 x lO^'^ for neutrino mass of ^ — 
1.8eV. This scale factor a„r is comparable with a^ec — 
9.2 X 10^** at recombination, but larger than acq ~ 10^^. 



D. Evolution of anisotropic stress 

From Eq. (pi]) the total anisotropic stress, tTj^^;' 
the r.h.s. of Eqs. and ((^ grows as iTtotai oc rf for 
both massless and massive neutrinos as long as the per- 
turbation is outside of the horizon and neutrinos are rel- 
ativistic. In the massless neutrino case, iTrl™^! decreases 
as time goes on and damps with some oscillations once a 
perturbation of wave number k enters the horizon. Then 



the PMF anisotropic stress 
at later epochs because t^'rmf 



dominates Itt. 



(m) 



'PMF\ uuiiiiiiauco 

remains constant. On 
the other hand, in massive neutrino case, the neutrino 



(^32) anisotropic stress |7r,'™''| decrease simultaneously at all 



scales of the perturbations if the wave number k satisfies 



k km, and |7r^,}^| quickly dominates {ttI'^^'^iI. For per- 
turbations at smaller scales with k ^ km, however, the 
neutrino mass effect is negligible and perturbations grow 
as if neutrinos have no mass. These features of massive 
neutrinos make large differences in the CMB. 



(m) 



where c = 45/14 for the vector mode (m — I), c — 
15/14 for the tensor mode (m = 2) and Hr is a density 
parameter of radiation when all neutrinos are relativistic. 
Inserting Eq. ([50)1 into Eq. pTjl . we obtain the initial 



E. PMF power spectrum 

Before closing this section we make a note on the power 
spectrum P{k) used in Eqs. (|30p and ([M)) . In the usual 
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CMB theory, the power spectrum is 

fc3p(fc) = A,fc"=-\ (36) 

where As is the scalar amphtude and Us is the scalar 
power spectral index. In the present calculations we need 
to use the power spectrum for the PMF [J, [3, i, i, i, ^ 
which can be written approximately as 

Pp(fc) oc /c2"B+6 . (37) 

It is to be noted that the accurate formula [l], [l^ of this 
rank of a single power spectrum has different forms for 
the scalar, vector and tensor modes for the PMF as a 
function of k and kc which is the cutoff wave number. 
See Ref. [l| and [l^l for more details. 



III. RESULTS AND DISCUSSION 

We have applied these equations which include a finite 
neutrino mass to the CMB anisotropy code CAMB 24] . 
As an illustration of the CMB anisotropics in the pres- 
ence of a PMF we consider an example with a fixed am- 
plitude and spectral index for the PMF power spectrum. 
Ultimately, one would hope to be able to deduce the am- 
plitude and spectral index from fits to the polarization 
spectrum. For our purposes, however, we assume that the 
PMF is generated before the nucleosynthesis epoch. In 
that case, the CMB anisotropics at high multipole places 
an upper limit on the magnetic field of Bx < 4.7nG [1|, 
and a nearly scale invariant spectrum is preferred [ll.[25|. 
Therefore, in the present work, we fixed the PMF at a 
field strength of B\ — 4.7nG, and adopt a spectral in- 
dex of riB = —2.9. As a representative neutrino mass, we 
choose a value of ^ rrii, = 1.8 eV which is the upper limit 
from the scalar mode analysis, '^m^ < 1.8 eV, deduced 
in [14], and is near the upper limit deduced from the 
WMAP-3yr analysis, "^"^ < 2 eV, 16]. Other param- 
eters are taken at the best fit values from the combined 
CMB, SNIa, LSS analysis from the WMAP-3yr data [H, 
except that the best fit value for flcDM is replaced with 
^CDM — r2^(1.8eV) in order to maintain a flat universe 
model. 



A. Power spectra 

Fig. [T] shows the power spectra for the scalar, vector 
and tensor, TT, EE, and TE modes along with the vector 
and tensor power spectra for the BB mode when a PMF 
is included. The upper most black lines, except for the 
regions 1^3 for TT mode and l^lOO for BB mode, show 
the primary spectra, and the other thin and thick lines 
represent models with massless and massive neutrinos, 
respectively, which include a PMF. 

Massive neutrinos have little effect on the scalar TT 
power spectra as evidenced by the fact that the thin 



green line and thick green line are almost indistinguish- 
able from each other. The only effect of massive neu- 
trinos in the scalar TT mode is a 3% enhancement for 
the quadrapole {I = 2) fluctuation. Also, the massive 
neutrinos make almost no difference in the vector (blue) 
and tensor (magenta) modes except at lower multipoles 
I^IOO. There, an excess of power at low I is caused by 
an increase in the shear as we now discuss. 



B. Analysis of shear 

To aid in this discussion. Figs. [5] and [3] show the ab- 
solute values of the shear a^™\ total anisotropic stress 
7r^™\ anisotropic scattering source P^™), and the ten- 
sor metric perturbation Hj^ . These quantities are nor- 
malized by the square root of power spectrum, P{k), 
as was defined below Eq. ([50]) . They are plotted as 
a function of the scale factor a for the vector (Fig. 
[5]) and tensor (Fig. [3]) modes. These quantities are 
shown for five different scales of the wave number k = 
10^^, 5 X 10-3, 10-3, 10-4, 10-5 Mpc-i from top to 
bottom in each panel of these figures. 

In each figure the l.h.s. panels display the models 
for massless neutrinos, and the r.h.s. panels display the 
models for massive neutrinos. The two upper-most lines 
for the larger wave numbers (smaller scales) {k = 10~^ 
and 5 x IQ-^ Mpc-^) are very similar to each other for 
both massless and massive neutrinos. However, the other 
larger scales k < lO-^Mpc-^ show a totally different evo- 
lution for massless vs. massive neutrinos. We understand 
the reason for similarity and difference of these quantities 
between the two models of massless and massive neutri- 
nos as we now describe. 



C. Growth of perturbations 

In the case when a perturbation is outside of the hori- 
zon, fcry <C 1, we know from Eq. ()34p that the normalized 
total anisotropic stress tt^^^; grows quadratically (oc rf). 
Using the fact that rj (x a during the radiation dominated 
epoch (a < aeq) and that rj cc a}!"^ during the matter 

dominated epoch (a > flgg), it is clear that tTj^^; oc 
(i.e. a = 2) in the radiation dominated epoch, and 

^'totai ^ ^ (^-^^ ck = 1) during the matter dominated 
epoch. This behavior is apparent in Figs. [2]and[3l 

In the same way we can understand the growth of 
shear. From Eqs. (P7)) and ([^5]) and a obtained from 
the above, the growth of the shear becomes a^™^ oc a 
during the radiation dominated epoch and 
in the matter dominated epoch, as is also evident in Figs, 
[landll 

As time goes on, the anisotropic stress of the neutri- 
nos decreases, and the PMF anisotropic stress becomes 
dominant, as can be seen in Eq. (|34p . However, once 
a perturbation enters the horizon, hierarchical mixing 
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sets in, and the initial condition Eq. (|34p is no longer 
valid. It is not easy to derive an analytic solution for the 
anisotropic stress. However, we know from numerical 
calculations that neutrino anisotropic stress undergoes 
damped oscillations, and tt^^q; remains constant asymp- 
totically as oc a° as shown in Figs. [2]and[3l Having this 
result, and applying a w to Eqs. (P7)) and we 
can expect (t*^™^ = cia"~^-^ + C2a^^ ~ cio"^^-^ oc a^^-^ 
during the matter dominated era. This damping power 
spectral index —1.5 is in good agreement with the slope 
of fco-(™) seen in Figs. [HandEl 

If the neutrinos have mass, the evolution of the 
shear ka^™-\ the total anisotropic stress tTj^^;, and the 
anisotropic scattering source p(™) drastically change as 
displayed on the r.h.s. panels of Figs. [2] and [3l In this 
case of finite mass neutrinos, initial conditions of the 
form of Eq. (j34p are valid when perturbations are out- 
side the horizon and when neutrinos behave relativisti- 
cally, i.e. when a < anr ~ 0.84 x 10""^. 

There is a distinctive feature of nearly equivalent evo- 
lution for these quantities especially at larger scales, 
fc<10~^Mpc"^ When one takes Y.'^v = 1.8eV, the 
critical wave number km of Eq. ([55)) is fcm 3.5 x 
10~'^Mpc~^ for the vector mode and k^ = 6.0 x 
10~'^Mpc"^ for the tensor mode. Therefore, k^ff ~ k 
for smaller scales, k > km, and k^ff « k„i — const. 
for larger scales, k < km- Putting these conditions into 
Eq. ([M)) , a similar evolution to that of massless neutrino 
models is expected for smaller scales k > km- However, 
the evolution becomes almost degenerate for larger scales 
k < km ^ 10~'^Mpc~^. This is the reason for the drastic 
change from massless neutrinos to massive neutrinos in 
ka^'^\ TT^^^^i and P^™) for larger scales with smaller wave 
number k < km- 

The effect of such changes brought on by massive neu- 
trinos is even more dramatic at low multipoles in the 
CMB power spectrum. The multipole Im corresponding 
to the scale km where the evolution becomes almost de- 
generate is 

Im kmVa , (38) 

where rjo ~ 14Gpc in the standard ACDM model. Since 
km is known to be 3.5 x lO-^Mpc"^ and 6.0 x 10~^Mpc"i 
for the vector and tensor modes, respectively, the critical 
multipole Im turns out to be 

Im ~ 50 (Vector mode) , 

Im - 85 (Tensor mode) . (39) 

Hence, the CMB power spectrum for lower mutipoles 
I < Im, corresponding to smaller k < fc,„, is expected 
to stay at the same value, which is in reasonable agree- 
ment with the calculated results shown in Fig. [1] In fact, 
there are also neutrino effects in the scalar mode at low I. 
However, they are ambiguous because of confusion from 
the integrated Sachs- Wolfe effect which is large when a 
PMF is present. 



D. EE and BB Modes 

The neutrino mass effect on the EE mode is larger 
than that of the BB mode. This difference is caused by 
the nature of the radial E and B functions el™"* {x) and 

Z?}™'' (x) . The behavior of these functions at smaller cc < 5 
is plotted in Fig. [4| for / = 2,3 and m = 1,2. From the 
definition of the radial E and B functions in Ref. [T9| . 
their leading order term behaves as ef^\x) oc x^~'^ and 
/3|™''(a;) oc x^^^ at x ~ 0, where we have used a form 
of spherical bessel function ji{x) oc at a; 0. Only 
the e"2^\x) fimction is finite at x = 0, i.e. e^'^ 0.2 
as X — > 0, while the other radial functions all vanish, 
— > at X ^ 0. This finiteness of 62™' (0), combined with 
a nearly scale invariant power spectrum for the PMF, 
i.e. k^Pik) oc fc^nij+e ^ ^0 j^j. ^ ^3 fj-om simple 
approximation Eq. (|37p . causes excess power for both 
the vector (m — 1) and tensor (m = 2) components in 
the EE and TE modes for / = 2, as we discuss below. 

Using the fact that the visibility function is approxi- 
mated by a delta function fe""^ ~ i5(ry — ijrec) and that 
OC fee// ^'^^ scales outside the horizon, we can esti- 
mate the CMB power spectrum of Eq. ([^0)1 as 

{21 + i)2cf 

« / ^k'^-^^Xfte^^'ikiVo-Vrec)) .(40) 

If neutrinos are light enough to be relativistic at recombi- 
nation, Eq. (|40p is a very good approximation for scales 
outside of the horizon for I'^IOO. Note that this approxi- 
mation is also valid in our calculation because anr ~ Orec- 
For the massless neutrino case, — 0, we can insert 
kf.ff = k into Eq. ([iU)) . Wc then obtain 

(2/ -f 1)2C72^^('") oc J dkk'--+h^r^'ik{r^o-Vr.c)) (41) 

This converges to a finite value for all multipoles I even 
when one takes a nearly scale invariant power spectrum 
P{k) with riB ^ —3. Our calculated CMB power spectra 
for the EE mode displayed in Fig. [1] does not show any 
excess at / < 10 in the massless case. However, if the 
neutrinos have mass, k^ff = y^k"^ + k'^ ~ km on very 
large scales, and Eq. ([iU]) becomes 

(2? + l)2cf^(™)(TO, ^0) 

oc J dkk^^-+^kie'i"'^\k{lJO'Vrec)) 

OC J (ifcfc2("«+3)+2'-5 ^ (42) 

where we have set e|'"''(x) oc x'~^ (x ~ 0) in the inte- 
grand on the r.h.s. When one takes a nearly scale in- 
variant power spectrum (ns ^ —3 ) for the PMF [P{k) 
of Eq. ([57)) ] , this integral has a logarithmic infrared di- 
vergence for the quadrapole I = 2 term, although it is 
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regular for the higher multipoles / > 3. In the present 
calculation, we have set = —2.9, because it is im- 
possible to set hb — —3 [1,11] although this nearly scale 
invariant spectrum is preferred P, [25| . This is the reason 
why our calculated EE mode shows a huge (but finite) 
excess for ^ = 2. 

For these mechanisms the anisotropics of the tensor EE 
mode becomes 100 times larger than the primary power 
spectrum. The ratio of massive to massless neutrinos, 
Ci^h/Ci^v, is a function of neutrino mass and spectral in- 
dex UB as we can see from Eq. (p4|) and ([42]) . The ratio 
does not depend on the amplitude of magnetic field, B\, 
because C; varies as C/ cx Q in the same manner in- 
dependently of the neutrino mass. Keeping this in mind, 
since strong enhancement of the anisotropics of tensor 
EE mode depends upon ub, and ^m^, it may be 
possible to observe this effect in the future and thereby 
place a strong constraint on all three of these quantities 
and on the neutrino mass in particular. 

There is however a large effect of cosmic variance in 
lower I, which is proportional to a/2/(2/ -|- 1) because of 
the finiteness of sampling [2^. This variance makes it 
difficult to clearly observe neutrino mass effect. Even if 
we cannot obtain lower limit of power spectrum due to 
the cosmic variance, we can still obtain the upper limit 
at lower I, from which we can constrain the upper limit 
of the neutrino mass. There is also a high probability of 
observing the neutrino mass effect in the BB mode if the 
gravitational wave is sufficiently weak, though we should 
consider carefully the passive mode studied by Lewis Q 
in BB mode. 
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Core-to-Core Program, International Research Network 
for Exotic Femto Systems (EFES). Work at UND sup- 
ported in part by the US Department of Energy under 
research grant DE-FG02-95-ER40934. 



IV. CONCLUSION 



In this work, we have expanded on earlier studies of 
the scalar CMB anisotropics in the presence of a PMF. 
In particular we have derived new vector and tensor mode 
equations in the presence of massive neutrinos and a 
PMF. We find a large effect from a finite neutrino mass on 
the vector and tensor modes when a PMF exists. In par- 
ticular, the effect of massive neutrinos on the EE mode 
become comparable to the observed primary anisotropy. 
Therefore, if and when the polarization power spectrum 
is ever measured at low multipoles, the possibility may 
exist to place a much stronger constraint on the sum of 
the neutrino masses than presently exists, though the 
effect of cosmic variance should be carefully taken into 
consideration. 
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FIG. 1: The power spectrum of the TT, EE, BB, and TE modes when a PMF is included. The black thin hues show the 
primary spectra, the scalar mode for the TT, EE, and TE mode, and the tensor mode for the BB mode. For this figure the 
amplitude ratio of the tensor-to-scalar mode has been set to 0.55. The other thin lines represent massless neutrino models. 
The thick lines show a model with massive neutrinos for which, rriv = 1.8 eV. Green lines are for the scalar mode, blue lines 
are for the vector mode, and red lines are for the tensor mode. The amplitude of the PMF is taken to be B\ = 4.7 nG, and 
the spectral index is taken to be ns = —2.9. Data points for the TT mode are from the WMAP-3yr data. 
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FIG. 2: Absolute values of the shear fca'^', total anisotropic stress ttHI^i, and the anisotropic scattering source P^^\ normalized 
by the square root of power spectrum, for the vector mode as a function of the scale factor a. The figures on the left hand side 
are for '^m,, = eV, and the right hand side figures are for ^^niy — 1.8 eV. Each line corresponds to a different scale, i.e. 
red: k = 0.01, green: k = 0.005, blue: k = 0.001, magenta: k = 0.0001, light blue: k = 0.00001 (Mpc~^). 
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FIG. 3: Absolute values of the tensor-mode metric perturbation h!^ , the shear ka^^\ the total anisotropic stress TftoLi' ^''^'^ 
the anisotropic scattering source P^-^\ normalized by the square root of power spectrum, as a function of the scale factor a. 
The figures on the left hand side are for ^in„ = eV, and the right hand side figures are for X] '^i' = 1-8 ^V. Each line 
corresponds to a different scale, i.e. red: k — 0.01, green; k = 0.005, blue: k = 0.001, magenta: k = 0.0001, light blue: 
k = 0.00001 (Mpc-^). 
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FIG. 4: Radial E function e['"'(a;) and radial B function /3['"'(a;), where Z = 2, 3 and m = 1,2. The left hand side is for the 
vector mode, and the right hand side is for the tensor mode. 
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